SOLVABILITY OF THE H"^ ALGEBRAIC RICCATI EQUATION IN 

BANACH ALGEBRAS 

AMOL SASANE 



O 

(N 



Abstract. Let _R be a commutative complex unital semisimple Banach algebra with the 
involution ■*. Sufficient conditions are given for the existence of a stabilizing solution to the 
H°° Riccati equation when the matricial data has entries from R. Applications to spatially 
distributed systems are discussed. 
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1. Introduction 

In the standard (regular full information) H°° control problem, one wants to find an appro- 
priate control input u such that the effect of the disturbance w on the output z is minimized 
for the control system given by the following dynamics: 

(1.1) x'{t) = Ax{t) + Bu{t) + Ew{t), 

(1.2) z{t) = Cx{t) + Diu{t) + D2w{t), 

where x{t) G C", u{t) G C™, w{t) G C^ z{t) € C^, and A,B,E,C,Di,D2 are complex 
matrices of appropriate sizes. More precisely, one seeks a feedback law 

u{t) = Fx{t) 

such that the closed loop system obtained in this manner is internally stable and its transfer 
function has H°° norm strictly less than some a priori given bound 7 > 0. The following 
result is well-known; see for example [HI Theorem 3.1, p. 49]: 

Theorem 1.1. Consider (jl.ip - (jl.2p . and let 'j > 0. Suppose that the tuple {A,B,C,Di) has 
no invariant zeros on the imaginary axis, and that kerDi = {0}. Then the following are 
equivalent: 

(1) There exists a matrix F such that the closed loop system obtained by applying the 
feedback u{t) = Fx{t) is internally stable, and the H^ norm of the closed loop transfer 
function is strictly less than 7. 

(2) D2D2 < 7^/. Moreover, there exists a positive semi-definite solution of the algebraic 
Riccati equation 

-1 



(1.3) = A*P+PA+C*C- 



B*P + DIC 
E*P + D^C 



DlDi 



DID2 



DID2 



I'l 



B*P + DIC 
E*P + D^C 



such that Aci is exponentially stable, that is, Re(A) < for all eigenvalues A of Ad, 
where 

(1.4) A,r.= A-[B E 



DtDi 



DID2 



D*2Di 



DID2 



Yi 



B*P + DIC 
E*P + D^C 



1991 Mathematics Subject Classification. Primary 46J05; Secondary 93D15, 58C15, 47N20. 
Key words and phrases. Riccati equations, Banach algebras, Systems over rings, H°° control. Spatially 
distributed dynamical systems. 



2 AMOL SASANE 

// P satisfies the conditions in (2) above, then the matrix F satisfying (1) can he taken as 

F ■= -{Dl{I - -i''^D2DlY^Diy^{DlC + B*P + DlD2{-f^I - DlD2r^{DlC + E*P)). 

The interest in systems over rings has recently been revived owing to potential applications 
for control theoretic problems of spatially invariant systems [1]. In particular, a natural 
question that arises is the following: if the matricial data in the algebraic Riccati equation 
has entries from a Banach algebra i?, then does there exist a solution matrix P which also has 
entries from Rl In the context of the algebraic Riccati equation associated with the problem 
of optimal control for a linear control system with a quadratic cost (LQ problem), this was 
done recently in [3j. In this article, we continue our study, and show that an analogous result 
can also be established for the H°° algebraic Riccati equation. Although the method of proof 
is identical to the one followed in [2] , the result in the current article does not follow from [3] , 
and it is not automatic. There are salient differences in the proof: 

(1) One of the key differences is that the quadratic term of the form P*QP in the algebraic 
Riccati equation has a positive semi-definite matrix Q in the LQ problem, whereas the 
matrix Q is indefinite in the H°° control problem. In particular, this manifests in that 
we need to prove the technical result (Proposition 12. 1( ) on the continuous dependence 
of the stabilizing solution of the H°° algebraic Riccati equation on the data. 

(2) The second new feature is the presence of the inverse of a matrix with feedthrough 
terms Di,D2 in the H°° algebraic Riccati equation. This results in an extra compli- 
cation in the application of the inverse function theorem for Banach algebras, namely, 
the application of Proposition 12.21 

We begin by fixing some notation. 

Notation 1.2. Throughout the article, R will denote a commutative, unital, complex, 
semisimple Banach algebra, which possesses an involution •*. 

On the other hand, the usual adjoint of a matrix M = [rriij] S £,px^ T^m \,q denoted by 
M* € C™^P, that is, M* = [m""]. 

M(R) will denote the maximal ideal space of R, equipped with the weak-* topology. For 
X G i?, we will denote its Gelfand transform by x, that is, 

x{(p) = ip{x), ip G M{R), X e R. 

For a matrix M G R^^"^, whose entry in the ith row and jth column is denoted by mij, we 
define M* G R'^^p to be the matrix whose entry in the ith row and jth column is m*,^. 

Also by M we mean the p x m matrix, whose entry in the ith row and jth. column is the 
continuous function m^j on M{R). Summarizing, if M = [iriij] G R^^"^, then 



mj. 



M* = 
M = 
[M{ip)j = ^i{^) 
Our main result is the following. 



G i?'"^^, 

pxm 



G [C{M{R);C) 

G O"""^. 



Theorem 1.3. Let A G i?"^", B G fl"^"*, C G iJ^^", Di G i?P^™, L>2 G PP""^ , E G i?"^^. 
Suppose that for all ip G M{R), 
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(Ai) iA^){^) = {A{ip)r, {B^){^) = {B{^)r, (c^)M = {c{^)r, {d{){^) = {D,{^)r, 
{dI){^1= {D2{v)yae^){^) = {E{v)r, 

(A2) {A{Lp),B{Lp),C{ip),Di{Lp)) is left invertible, 

(A3) there exist matrices Fi{ip),F2{(p) such that A{(p) + B{ip)Fi{ip) is exponentially stable, 
and we have (suppressing the argument ip) 

\\{C + DiFi){sI -A- BFi)-\E + BF2) + D2 + DiF2\\oo < 7, 

(A5) {A{ip),B{ip),C{{p),Di{ip)) has no invariant zeros on the imaginary axis, and 
(A6) ker(^i((^)) = {0}. 
Then there exists a P ^ j^nxn g^^/^ ^/jq,^ 

(1) = A*P + PA + C*C 

-[ PB + C*Di PE + C*D2 

(2) The matrix 



D\Di D\D2 



1 -1 



B*P + D*^C 
E*P + D^C 



D\Di 



D\D2 



D^Di D^D2 - YI 



-1 



B*P + D\C 
E*P + D^C 



A,i:=A-[B E] 

is exponentially stable, and 
(3) for all (f G M{R), P{(p) is positive semidefinite. 



In the following we define what is meant by "exponentially stable" in the conclusion (2) in 
the above result. 

Definition 1.4. Let i? be a commutative, unital, complex, semisimple Banach algebra. If 
A € i?"^", then let M^ '■ R^ — >■ R^ be the multiplication map by the matrix A, that is, 
V >—^ Av {v G R^). Then R^^^ is a unital complex Banach algebra (for example) with the 
norm 

\\A\\:=\\MA\\ciR-) (yl€ii"X") 

where C{R^) denotes the set of all continuous linear transformations from i?" to i?", and R^ 
is the Banach space equipped (for example) with the norm 



maxlllxfcll : 1 < A; < n} for X 



Xl 



and C{R^) is equipped with the usual operator norm: 

\\Ma\\c(R^) = suplll^ull : t; G i?" with \\v\\ < 1}. 

For A G R^^"^, we define 

00 ^ 



fc=0 



The absolute convergence of this series is established just as in the scalar case. 

A G R^^"- is said to be exponentially stable if there exist positive constants C and e such 
that 



JA\ 



< Ce"'* for ah t > 0. 
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We recall the following two results [U Lemma 1.6, Prop. 1.7]. 

Lemma 1.5. Let A € i?"^". Then the following are equivalent: 

(1) A belongs to the spectrum of A (^ i?"^"". 

(2) A belongs to the spectrum of Ma G C{R^). 

(3) A belongs to the spectrum of A{ip) for some ip € M{R). 

Proposition 1.6. Let A G i?"^". Then A is exponentially stable if and only if 
sup{Re(A) : X is an eigenvalue of A{{p) for some ip € M{K)\ < 0. 

The proof in Section [2] of our main result above is established by taking the Gelfand 
transform of our equation, and showing that the pointwise solution is continuous. Then we 
use the Banach algebra operational calculus to ensure that this continuous solution is actually 
the Gelfand transform of a matrix with entries from the Banach algebra. 

In Section [3] we discuss the applications of this result to the control of spatially invariant 
systems. 

2. Proof of the main result 

We will need the following two results. The first one gives sufficient conditions on the 
matricial data for the classical H°^ algebraic Riccati equation to have a positive semi-definite 
stabilizing solution, and moreover it says the smallest such solution P depends continuously 
on the matricial data. 

Proposition 2.1. Let S be the set of matrix tuples {A,B,C,Di,D2,E) belonging to C"^** x 

(^nxm ^ i^pxn ^ i^pxm ^ i^pxl ^ i^nx£ ^^^f^ ^^^^ 

(1) {A,B,C,Di) is left invertible, 

(2) there exist matrices Fi,F2 such that A + BFi is exponentially stable, and 

||(C + DFi){sI -A- BFiY\E + BF2) + D2 + ^iFalU < 7, 

(3) {A,B,C,Di) has no invariant zeros on the imaginary axis, and 

(4) kerL>i = {0}. 

Then there exists a smallest positive semidefinite, solution P = P{A, B, C, Di,D2, E) G (^"X" 
of the Riccati equation (jl.3p such that A^i given by (jl.4p is exponentially stable. Moreover, 
the map {A,B,C,Di,D2, E) ^ P{A, B,C, Di, D2,E) : S ^ C"^" is continuous. 

Proof. The proof is the same as that of [H Lemma 3.1], mutatis mutandis, except we have 
complex matrices instead of real ones, and transposes there have to be replaced by the Hermit- 
ian adjoints. We give repeat the proof here. The existence and uniqueness of the stabilizing 
solution is well-known; see [8]. The continuity follows in a straightforward manner. The 
solution P is associated with the stable subspace of a Hamiltonian matrix . Since this Hamil- 
tonian matrix has no eigenvalues on the imaginary axis, the eigenvalues in the open left half 
plane and the open right half plane are separated, and the existence of a continuous basis for 
the stable subspace and hence the continuous dependence of the stabilizing solution of the 
H°° algebraic Riccati equation can be found for example in [7]. D 

The next result we will need is the following version of the Implicit Function Theorem in 
the context of Banach algebras (see [H p. 155]), and this will be used in proving our main 
results when we need to pass from continuous functions on M(R) to elements of R. 
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Proposition 2.2. Let hi, . . . ,hs be continuous functions on M{R). Suppose that fi, . . . fr in 
R and Gi{zi, . . . , z^^r), ■ ■ ■ , Gt{zi, . . . , Zg+r) <ire holomorphic functions with t > s defined on 
a neighbourhood of the joint spectrum 



a{hi,... ,hs,fi,...,fr) 
such that 
(2.1) Gk{hu...,hs,fi,...,fr) 

If the rank of the Jacobi matrix -— — ^ ^ — — 

0{zi,...,Zs) 

elements gi, . . . ,gs in R such that 



{(hiiip), ..., h,{^), fiiip), ..., fri^)) : ^ G MiR)}, 

^ =0 on M{R) forl<k <t. 

is s on a{hi, . . . ,hs, fi, ■ ■ ■ , fr), then there exist 



gi =hi,...,gs 
We are now ready to prove our main result. 



h... 



Proof of Theorem \1.3[ If we fix a (/? G M{R), then owing to the assumptions (A2)-(A6), we 
know that there is a unique smallest positive semidefinite solution in C"^", which we will 
denote by 11(f), such that (suppressing the argument if) 



(2.2)0 = {A)*n + UA + {C)*C 

{B)*U + {D~i)*C 

{E)*U + {D2)*C 



{Di)*Di {Di)*D2 

(D2)*Di (D2)*D2-^^I 



{B)*U + {Di)*C 
{E)*U + {D2)*C 



and 



A 



cl 



A-[ B E 



iDi)*Di {D1YD2 

(D2)*Di {D2)*D2-^H 



{Byn + iDiYc 

{E)*\l + {D2)*C 



is (pointwise, for each tp G M[R)) exponentially stable. 

Moreover, from Proposition 12. 11 it follows that the map ip 1-^ 11(99) is continuous on M{R). 
In light of the assumption (Al) and (12. 2p . we can hence conclude that there is a continuous 
solution n on the maximal ideal space such that 







{A'')U + UA + {C*)C 

UB + C*Di UE + C*D2 



D\Di D\D2 

dIDi D*D2-j^I 



B^'U + D^C 
E'*U + D*d 



Finally we will apply Proposition 12.21 We have in our case s = n^, t = r? , the /ij:s are the 
components of 11 and the /j:s are the components of A, A*' , B, B*, C, C* , Di, D\, D2, D2, E, E* 
(which are totally r = 2(n^ + nm + pn + pni-\-pl + ntj in number). The maps Gi,. . . , G^=„2 
are the n^ components of the map Q 

{e,Ui,U2,Vi,V2,Wi,W2,Xi,X2,Yi,Y2,Zi,Z2) 



[ QVi + W2X1 ezi + W2Y1 



U2Q + QUi + W2W1 

X2A'i X2Y1 

Y2XI Y2YI - 72/ 



V2Q + X2WI 
Z2G + Y2WI 
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(In the above, we have the replacements of 

A,A\B,B\C,C\Di,D{,D2,D*,E,E* 
by the complex variables which are the components of 

Ui,U2,Vi,V2,Wi,W2,XuX2,Yi,Y2,Zi,Z2, 

respectively. The replacements of the components of 11 in the Riccati equation are by the 
complex variables which are the components of 0.) Since the set of invertible matrices is 
open, it follows that 

{X,,X2,Yi,Y2)^ 



X2X1 X2Y1 

Y2X1 Y2Y1 - 72/ 



is holomorphic on the joint spectrum of {Di,D*,D2,D2)- Consequently the components 
Gi, . . . , G^2 of above map Q is holomorphic on the joint spectrum. 

There is a continuous solution 11 on the maximal ideal space such that for all A;, 



(2.3) Gfc(n, A, A-, B, B\ C, C*, D^,D\, D2, D*, E, E*) = 

on M{R) (that is, condition (j2.ip in Proposition 12.21 is satisfied). So we now investigate the 
Jacobian with respect to the variables in Q. The Jacobian with respect to the variables at 
the point 

(n((^), I(^), {A{^))\ B{^), {B{^)r, c{ip), {d{^)r, 

2 2 

can be verified to be the following linear transformation A from C" — )■ C"" : 

The set of eigenvalues of A consists of the numbers 

-(A + /i), 

where A, fi belong to the set of eigenvalues of Aci{^); see for example [21 Proposition 7.2.3]. 
But since Aci{f) is exponentially stable, all its eigenvalues have a negative real part. Hence 
— (A+/Li) 7^ for all A, fj, belonging to the set of eigenvalues of the matrix Aci{^)- Consequently, 
the map A is invertible from C" to C" , and its rank is n? = s. So by Proposition 12.21 there 
exists a P G i?"""" such that P{ip) = n(c^) for all if £ M{R). From ([23]), it follows (using the 
fact that R is semisimple) that 

= A*P + PA + C*C 

D\Di D\D2 1 "^ r B*P + D\C 

DlDi DID2 - 7^/ J [ E*P + L>JC 

From the property possessed by the pointwise solutions n((/7) {ip £ M{R)) of the constant 
complex matricial Riccati equations (j2.2p . we have that for all 99 G M{R), all eigenvalues of 
Acii^p) have a negative real part. But the set-valued map taking a square complex matrix of 
size n X n to its spectrum (a set of n complex numbers) is a continuous map; see for example 
[6l II, §5, Theorem 5.14, p. 118]. Since M[R) is compact in the Gelfand topology (the weak-* 
topology induced on M{R) considered as a subset of C{R; C)), it follows that 

sup{Re(A) : A is an eigenvalue of Ac\{<^) for some ip G M[R)} < 0. 



[ PB + CDi PE + C*D2 
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From Proposition [121 it follows that A^y is exponentially stable. Finally, again by the property 
possessed by the pointwise solution 11, we have that for all ip € M{R), P{ip) is positive 
semidefinite. This completes the proof of Theorem 11.31 D 

We observe that whether or not the assumption (Al) in Theorem 11.31 holds is intimately 
related to the choice of the involution •* in the Banach algebra R. For some commutative 
Banach algebras with involutions, this is automatic, namely if it is symmetric; see [4]. 

3. Application to spatially invariant systems 

Finally, we mention the application of our main result to control problems for spatially 
invariant systems introduced in [1]. The analysis of spatially invariant systems can be greatly 
simplified by taking Fourier transforms, see [1], [3]. This yields systems described by mul- 
tiplication operators with symbols belonging to L°°(T). The H°° control design is to use 
the feedback F, which uses the bounded, self-adjoint, stabilizing solution P to the H°° al- 
gebraic Riccati equation on the Hilbert space (L^(T))". For the design of implementable 
controllers it is important that the gain operator have a spatially decaying property (see [1]). 
This translates into the mathematical question of when the H°° algebraic Riccati equation 
has a stabilizing solution in a suitable subalgebra (for example, L^(T) is a subalgebra of 
L°°{T) C C{L'^{T))). So the spatially decaying property now translates into finding suitable 
subalgebras of L°°{T), in particular, the weighted Wiener algebras. In [H p. 472-474], an 
example of a large class of the Wiener subalgebra of L°°{T) having the symmetry property 
(which leads to automatic satisfaction of the condition (Al) in our main result) was given. 

Acknowledgements: The author thanks Ruth Curtain for suggesting to the author the 
question of investigating if the result from [4] is valid also for the H°° algebraic Riccati 
equation, and also for useful discussions on the classical H°° problem. 
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